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We describe a scalable stochastic method for the experimental measurement of generalized fidelities charac- 
terizing the accuracy of the implementation of a coherent quantum transformation. The method is based on the 
motion reversal of random unitary operators. In the simplest case our method enables direct estimation of the 
average gate fidelity. The more general fidelities are characterized by a universal exponential rate of fidelity 
loss. In all cases the measurable fidelity decrease is directly related to the strength of the noise affecting the 
implementation - quantified by the trace of the superoperator describing the non-unitary dynamics. While the 
scalability of our stochastic protocol makes it most relevant in large Hilbert spaces (when quantum process to- 
mography is infeasible), our method should be immediately useful for evaluating the degree of control that is 
achievable in any prototype quantum processing device. By varying over different experimental arrangements 
and error-correction strategies additional information about the noise can be determined. 

INTRODUCTION 

The characterization and elimination of decoherence and other noise sources has emerged as one of the major challenges 
confronting the coherent experimental control of increasingly large multi-body quantum systems. Decoherence arising from 
undesired interactions with background (or environment) systems and imprecision in the classical control fields lead to severe 
limits on the observation of mesoscopic and macroscopic quantum phenomena, such as interference effects, and, in particular, the 
realization of quantum communication and computation algorithms. Measurement of the strength and other detailed properties 
of the noise mechanisms affecting a physical implementation is a critical part of optimizing, improving, and benchmarking 
the physical device and experimental protocol ma. Moreover, in the case of quantum devices capable of universal control, 
knowledge of specific characteristics of the noise enables the selection and optimization of passive and active error-prevention 
strategies 013 00 011]. 

The exact method for characterizing the noise affecting an implementation is quantum process tomography (QPT) 1 8]. Let D 
denote the dimensionality of the Hilbert space (constituted, e.g., from n q = log 2 (D) qubits). For QPT, the desired transformation 
(usually a unitary operator) must be applied to each member of a complete set of D 2 input states (spanning the state space), 
followed by tomographic measurement of the output state. This allows for a complete reconstruction of the superoperator 
(completely positive linear map) representing the imperfect implementation of the target transformation. From this superoperator 
the cumulative noise superoperator can be extracted from conventional analysis of the matrix. The QPT approach to noise 
estimation suffers from several practical deficiències. First, often the intrinsic properties of the noise operators are of interest, 
but the noise superoperator determined from QPT will depend on the symmetries between the noise mechanisms and the choice 
of target transformation. Second, the number of experiments that must be carried out grows exponentially in the number of 
qubits D 4 = 2 4rlq . Third, conventional numerical analysis of the tomographic data requires the manipulation of matrices of 
exponentially increasing dimension (D 2 x D 2 ). For these last two reasons QPT becomes infeasible for processes involving more 
than about a dozen qubits, far fewer than the one thousand or so qubits required for the fault-tolerant implementation of quantum 
algorithms that outperform conventional computation. Hence the infeasibility of complete noise estimation via tomography 
prompts the question of whether there exist efficient methods by which specific features of the noise may be determined. 

We show below that the overall noise strength and the associated accuracy of an implementation may be estimated by a 
scalable experimental method. Specifically we show that the average gate fidelity Q, and some more generalized fidelities 
described below, can be estimated directly with an accuracy 0(l/yDN) where N is the number of independent experiments. 
This method provides a solution to the important problem of efficiently measuring which member of a set of experimental 
configurations and algorithmic techniques produces the most accurate implementation of an arbitrary target transformation. By 
varying over different experimental methods and noise-reduction algorithms and then directly measuring the variation in the 
associated fidelity this method enables estimation of more detailed characteristics of the noise. 
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EFFICIENT ESTIMATION OF THE AVERAGE GATE FIDELITY 



A convenient starting point for our analysis is the average gate fidelity 

F^(A) = E^(F e (U,A,^))= /<# (W~HMU\Í>)MU- l )m) (D 

where 

A(p)=J2 A kp4 (2) 

fc 

is a completely positive (CP) map characterizing the noise. The gate fidelity F g is the inner-product of the state obtained from 
the actual implementation with the state that would be ideally obtained under the target unitary. The measure dip denotes the 
natural, unitarily invariant (Fubini-Study) measure on the set of pure states and hence the average gate fidelity provides an 
indicator that is independent of the choice of initial state. If the implementation is perfect then F g = 1 and under increasing 
noise F g decreases. Due to the invariance of the Fubini-Study measure the average fidelity depends only on the noise operator 
and can been expressed in the form I9t ll0lll 111 

Hence the average fidelity can be determined if the noise operator is known. The noise operator can be determined experimentally 
by measuring the CP map A{U ■ U^ 1 ) tomographically and then factoring out the inverse of the target map U^ 1 ■ U. This 
procedure has been carried out recently for 3 qubits in recent a implementation of the quantum Fourier transform using liquid- 
state NMR techniques |2J]. As noted above, this method requires 0(D 4 ) experiments and the conventional manipulation of 
matrices of dimension D 2 x D 2 . Recently Nielsen has proposed a method 1 11] for the direct measurement of F g that requires 
D 4 experiments but analysis of matrices of dimension only D x D (rather than D 2 x D 2 ). 

We now describe how the average gate fidelity {3 can be estimated accurately from a simple experimental protocol. Our 
method requires the physical implementation of the "motion reversal" transformation U~ 1 U\i(>) (ip\U~ 1 U on an arbitrary state 
|'í/>)(i/'|· Under this transformation, the CP map A in the gate fidelity Q can be interpreted as the decoherence and experimental 
control errors arising under the imperfect implementation of the motion reversal experiment, i.e., A = Ay-iy, rather than as 
the noise associated with only the forward transformation U, i.e, A = Ajj. The key idea is to choose the target transformation 
U randomly from the Haar measure 1 12]. This will earn us the advantage of the concentration of measure in large Hilbert 
spaces, as described further below, and leads to a universal form of the gate fidelity depending only on the intrinsic strength of 
the cumulative noise. This universal form will allow us to evaluate the average fidelity for more generalized motion reversal 
protocols. 

Our starting point is the gate fidelity uniformly averaged over all unitàries, 

^u{F g )= í d[/Tr[pC/- 1 A(í/pf/- 1 )í7], (4) 

J U(D) 

where in the above dU denotes the unitarily-invariant Haar measure on U(D) and p = \ip) (tp\. In order to evaluate this integral 
we use the superoperator representation of the map (|2jl, 



A = J2 A k®A* k , (5) 



fc 

and similarly Ü — U ® U*, where * denotes complex conjugation. The Haar averaged gate fidelity takes the form 

= i r | f> i dU UAU~ 



Pj (6) 
Tr (M avc p) = ^g(A avc ). (7) 



where À avc = J dU ÜAÜ 1 . As shown in the Appendix, the Haar-averaged superoperator À avc is [/(-D)-invariant and thus can 
be expressed as a depolarizing channel 

À a ™p = w +(l-p)l (8) 
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(assuming Tr(p) = 1) characterized by the single "strength" parameter 

= EJTr(A fc )| 2 -l 

P D 2_l 

where p e [0, 1] and we have made use of the fact that Tr(A avc ) = Tr(A) = J2k I^O^fc)! 2 - Direct substitution leads to 



(9) 



Eu(F s ) = F g (À avc ) = p + íi^. (10) 

Hence the gate fidelity for the Haar-averaged operator resulting from a motion reversal experiment depends only on the single 
parameter Tr(A) which represents the intrinsic strength of the cumulative noise. We remark that this result holds for general 
(possibly non-unital) noise. Furthermore, suppressing the arguments of F we note that the unitary invariance of the natural 
measure on pure states implies the equivalence 

F r g (A)=E i ,(F g )=E u (F g ), (11) 

and hence we recover Eq.|3] 

We now describe why and how the intrinsic noise strength (characterized by p or Tr(A)) can be estimated via an efficient 
experimental protocol. By implementing a single target transformation U that is randomly drawn from the Haar measure, we 
gain the advantage of the concentration of measure in large Hilbert spaces: the motion reversal (gate) fidelity for the single 
random U is exponentially close to the Haar-averaged motion reversal (gate) fidelity. From the unitary invariance of the Fubini- 
Study measure we know that 

E^F*) = E U (F*). (12) 

As will be shown in Ref. 1 13], the typical fluctuation for a random initial state \ip), given a fixed U and A, decreases exponentially 
with the number of qubits, 



(A^) 2 = - F g 2 < 0(1/D). (13) 



Therefore it follows that, 



(AFYu = Eu(F') - Eu(F g Y < 0(1/ D). (14) 



Hence the fidelity under motion reversal of a single random U and arbitrary (non-random) initial state is exponentially close to 
the Haar-averaged fidelity 

F S (U, A, = F g (A ave ) + 0(1/ VD) =P+ + 0(1/VD). (15) 

The protocol is now clear: after the motion reversal sequence has been applied experimentally, the single parameter p charac- 
terizing the average gate fidelity appears as the residual population of the initial state. Due to the invariance of the Haar measure 
we may choose the initial state to be the computational basis state (|0)(0|)® n<! . Hence the gate fidelity can be determined di- 
rectly from a Standard readout (projective measurement) of the final state in the computational basis. When the noise strength 
is actually non-negligible (e.g., the noise strength does not decrease as a polynomial function of 1/D) an accurate estimate of 
p is possible with only a few experimental trials. If in each of N repetitions of the motion-reversal experiment an independent 
random unitary is applied, then the observed average will approach the Haar-average as 0(l/yDN). 



GENERALIZED FIDELITIES IN A DISCRETE-TIME SCENARIO 



More generally we imagine the ability to implement a set of independent random unitary operators {Uj} and their inverses. 
The entire sequence is subject to some unknown noise, consisting of the decoherence processes and control errors affecting the 
implementation. Such generalized motion reversal sequences are relevant not only for noise-estimation, but also have important 
applications in studies of fidelity decay 1 14] and decoherence rates for quantum chaos and many-body complex systems. 

We first consider the fidelity loss arising under an iterated motion reversal sequence of the form 



p(n) = Ü^kÜn . . . Ü^kÜ2 Ú^AÜ lP (0), 



(16) 
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where here Aj = À^-i^ denotes the cumulative noise from the motion reversal of Uj and we now allow arbitrary (possibly 
mixed) initial states p(0). The fidelity of this iterated transformation is, 



F„(V>, {Uj}) = Tr (p(0)[/ I 7 1 À„í7„ . . . Ü^ktÜípioj) . (17) 
Averaging over the Haar measure for each Uj takes the form, 

ï^ = IE {a . } (F„(V>, {£/;})) = í (Tl% x dUj)F n {^{Uj}) (18) 



Trfp(O) H? =1 Àf e /,(())). .19) 



where dUj denotes the Haar measure and we have defined the Haar averaged noise operator, 



Aj vc = JEjjj (Aj) = / dUU- l kjU. (20) 

JU(D) 



As noted above and shown in the Appendix, À avc = j dUÜ AÜ 1 is a depolarizing channel 



Af> = ft p+(l- Pj )I (21) 

with strength parameter 

Because each Uj is random, we can further simplify this result by assuming that the cumulative noise for each Uj has the same 
strength pj = p, in which case we obtain for arbitrary noise a universal exponential decay of the averaged fidelity 

7^ = p"Tr[p(0) 2 ] + ^-^. (23) 

depending only on the noise str eng th. In the limit of large n, we see that F n — ► £> _1 , as may be expected from the average 
fidelity between random states 11611 . Most importantly, due to the concentration of measure J14-I - for large D the fidelity loss 
under iterated motion reversal of a single sequence of random unitary operators will be exponentially close to the Haar-average, 
and hence the noise strength can be estimated with only a few experimental runs. 

Another important generalized fidelity is the one obtained under the imperfect 'Loschmidt echo' sequence 

p(n) = Ü^ 1 . . . Ü^AjJn . . . Ài#ip(0), (24) 

where the superoperator Aj represents the cumulative noise during the implementation of each Uj. The fidelity between the 
initial state and final state in the Loschmidt echo experiment takes the form, 

í; ccho 0A, {A,}, {Uj}) = Tr {p{Q)Ü^ l Ü^ . . . Ü-^Ajïn . . . 1^(0)) . (25) 

Moving to the interaction picture we define 

Aj{j) = Üï 1 ...Ür l AjÜj...Ü 1 , (26) 

so that, 

í; ccho (V, {Aj},{Uj}) - Tr (p(0)A n (n)A n _!(n - 1) . . . Ai(l)p(0)) . (27) 
From the invariance of the Haar measure the average fidelity simplifies to 



F echo = Tr /9 (o)A^ e A^_ e 1 . . . Af c p(0) (28) 
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with Aj vc given by Eq. |^ As before, we can simplify this result by assuming that the cumulative noise for each step has the 
same strength (pj — p), in which case we obtain for arbitrary noise a universal exponential form for the decay of fidelity 



jy ho (p)=p n + V £ ' ■ (29) 

A generalized version of this Loschmidt echo that is more relevant to noise estimation is one for which noise appears in both 
the forward and backward sequence of the motion reversal. The associated fidelity is, 

Ff A, {Ui}) = Tr (p{Q)kÜ^kÜ^ . . . kÜ^kÜn . . . Àüip(0)) . (30) 

While we have not directly evaluated the average of this fidelity analytically in the general case, for the special case of unitary 
noise we have analytic and numèric al evidence supporting the relation 

F gen „ ^ccho (31) 

for large n, which we conjecture should hold under general noise. 



GENERALIZED FIDELITIES FOR CONTINUOUS-TIME WEAK NOISE 

We describe our system by the Markovian Master Equation 

j t P = -i[Hc{t),p]+eL{p) (32) 
where Hc(t) governs a controlled reversible part of the dynamics and the generator 

L p = L(p) = -i[H, p] + ^([V a , pVÍ\ + [V a p, V^}) (33) 

a 

with the condition Trií = TrV a = (which fixes the decomposition of L into Hamiltonian and dissipative parts [ 19]) describes 
all sources of imperfections and noise. Here < e <C 1 is a small parameter characterizing noise strength. 
The time dependent fidelity of the initial state <f> is given by 



F^t) = (0|Texp{e J L(s)ds] (\cf>) (cl>\)\<t>) (34) 
where T denotes the chronological order, and 

L(s) = Ü Jt (s,0)LÜ{s,0), U(t,s)=Texp^-iJ Hc(u)duj. (35) 

Using the notation 

f(í) = Texpje^ L(s)ds} (36) 
we can write down the following "cumulant expansion" of the dynamics with respect to the small parameter e 

f(t) = exp{e^!(í) + e 2 K 2 {t) + •••}. (37) 
Using the Wilcox formula for the matrix-valued functions 



dx 

one obtains 



^ exp A(x) = i^J exp(\A(x))-^-A(x) ex_p(~\A(x))d\j exp A(x) (38) 



1 '* 



K l {t) = I L(s)ds, K 2 (t) = ~l ds I du[L(s), L(u)} , (39) 
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We assume now the following ergodic hypothesis: a) the ergodic mean exists and is equal to the Haar average 



1 ' T 



JU(D 



Hm - / L(t)dt = L avo = / dU ULU^ , (40) 



b) the fluctuations 8L(t) = L(t) — L avc around ergodic mean are normal, i.e. for long t 

rs+t 



(42) 



SL(u)du\\ ~ í 1/2 . (41) 

These condi tions are satisfied , for instance if the time-dependent dynamics 1 1— > Í7(í) can be modelled by a random walk on the 
group U(D) or by a trajectory on U(D) given by a certain deterministic dynamics with strong enough ergodic properties. The 
norm of ^(í) can be estimated using Í411 

\\K 2 (t)\\ = i || f ds J* du([5L(s),5L(u)] + [5L(s),L ave ] + [L^,SL(u)})\\ ~ í 3 / 2 . 
Therefore for small enough e and long enough times í such that ei is fixed the first term dominates and we can write 

f (í) ~ exp(e [ L(s)ds) . (43) 
Jo 

Then replacing T(t) by exp(eL av t) and using the explicit expression A56I57I we obtain the universal exponential decay of the 
fidelity 

F,(t),^ + i(l-^), 7= 2(J f_l) E tr (i^l 2 )- ^ 
DISCUSSION 

We have described how generalized Haar-averaged fidelities may be directly estimated with only a few experimental mea- 
surements. By implementing a motion reversal sequence with a Haar-random unitary transformation, the observed fidelity decay 
provides a direct experimental estimate of the intrinsic strength of the noise. Moreover, because the target transformation is a 
Haar-random unitary, the cumulative noise measured by this method will not be biased by any special symmetries of the target 
transformation. 

The only inefficiency of our protocol is the requirement of experimentally implementing a Haar-random unitary: the decom- 
position into elementary one and two qubit gates requires an exponentially long gate sequence |12J. However, the randomization 
provided by Haar-random unitary operators may be unnecessarily st rong an d this leads to the open question of whether efficient 
sets of random unitàries, e.g. the random circuits studied in Refs. J20l.l2 lli , can provide an adequate degree of randomization 
for the above protocols. Indeed the experimental results of Ref. |2] suggest that even a structured transformation such as the 
quantum Fourier transform is sufficiently complex to approximately average the cumulative noise to an effective depolarizing 
channel, and from studies of quantum chaos it is known that efficient chaotic quantum maps are faithful to the universal Haar- 
averaged fidelity decay under imperfect motion-reversal 0- While more conclusive evidence is needed to answer this question, 
it appears likely that the inefficiency associated with implementing Haar-random unitary unitary operators may be overcome. 

An additional question is whether the implementation of random unitary operators (e.g., Haar-random unitary operators or 
even efficient random circuits) leads to an even stronger form of averaging. We have throughout our analysis made the usual 
assumption that the noise superoperator A is independent of the specific target transformation but depends only on the duration 
of the experiment. However it is known that the actual noise in general depends sensitively on the choice of target transformation 
U. Moreover, the cumulative noise operator generally also depends on the particular sequence of elementary one and two qubit 
gates applied to generate U. For example, the implementation of the quantum Fourier transform lli |2| will generate very 
different cumulative noise than the trivial implementation of the identity operator U — 1 for the same time r. However, it 
appears likely that the cumulative noise operators, and in particular their intrinsic noise strength, under a specific but random 
gate sequence should become concentrated about an average value depending only on the length of the sequence. If this is the 
case, then the usual assumption that the noise is independent of the actual gate sequence becomes statistically well motivated, 
and the measured fidelity under motion reversal can provide a benchmark of an intrinsic noise strength that is fully independent 
of the target unitary. 

Note added in proof: additional evidence for the conjectured relation (31) can be found in Ref. j22|. 
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APPENDIX: HAAR AVERAGED SUPEROPERATORS 



We consider a linear superoperator À acting on the space Mjj of D x D complex matrices treated as a Hilbert space with a 
scalar product (X, Y) = Tr(X^Y). The superoperator A has a D 2 x D 2 dimensional matrix representation and TrA denotes 
the usual sum over the diagonal elements of the matrix. For clarity of notation we will sometimes express the linear operation 
Ap in the form A( ( o). By {|fc)} we denote an orthonormal basis in C D while {Eki = is a corresponding basis in M^». 

The group U(D) of unitary D x D matrices has its natural unitary representation onM^ defined by 

U{D) 3 U i ► U , ÜX = UXU f . (45) 

This representation is reducible and implies the decomposition of into two irreducible invariant subspaces 

M D =M C D ®M° D , M.% = {X G M-d; TrX = 0} , M C D = {X = cl}, (46) 



where c is an arbitrary complex number. 

Any superoperator A possesses exactly two linear U (D) invariants, i.e. the linear functionals on superoperator space which 
are invariant with respect to all transformation of the form A i— > UAU^ : 

D 

Tr[A(l)] = ^(fc|A(l)|fc) (47) 
k=l 

and 

D D 

Tr(À) = (Ekl,A(Ekl)) = ( k \MEki) \l) (48) 

k,l=l k,l=l 

Example Take A(X) = AXB, then Tr[A(l)] = Tr(AB) and Tr(À) = Tt(A)Tt(B). 

A U(D) -invariant operator satisfies A mv = UA^W for any U S U(D). The following lemma completely characterizes 
U( 7J)-invariant trace-preserving superoperators 

Lemma 1 Let A lnv be a U (I?)-invariant trace-preserving operator. Then 

À inv X = A ïm {X) =pX + (l-p) Tr(X)-^ , (49) 

where 

Tr(À inv ) - 1 

P = D 2 -l ■ (50) 

Proof Schur's lemma implies the form ( I49l > for Í7(l?)-invariant trace-preserving operators. From the normalization 
Tr[A lnv (l)] = D for the trace, the detailed expression (1501 can be explicitly calculated by comparing í7(D)-invariants for 
both sides of eq.(|49j. □ 

The Haar-averaged superoperator corresponding to the noise under the imperfect motion-reversal protocol, averaged over all 
possible unitary operators, is a U (£>)-invariant superoperator 

À avc = / dU ÜAÜï . (51) 

Ju(D) 

where dU is the normalized Haar measure on U(D). Using Lemma 1 we can easily compute the averaged form of the dynamical 
map for both the Schròdinger operator 

A(p)=Y,A aP Al, 5>tA Q = l (52) 

and for the semigroup generator 

L p = L(p) = -i[H, p] + pVl] + [V aP , V^]) (53) 
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with the condition Trií = TrV a — which fixes the decomposition of L into Hamiltonian and dissipative parts. From the fact 
that Tr(À avo ) = Tr(À) we obtain 

A avc (p)=pp+(l-p)Tr(p)l (54) 



where 



Similarly for the generator we obtain 



where 



Tr(A)-l = EJTr(A a )| 2 -l 
P D 2 -l D 2 -l ' K ' 



L'"' r ^L-U>) = - : [p-Tr(p)^) (56! 



^TyE^I 2 )- (57) 
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